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Radiative improvement of the lattice NRQCD action using the background field method and 
application to the hyperfine splitting of quarkonium states 
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We present the first application of the background field method to Non-Relativistic QCD (NRQCD) 
on the lattice in order to determine the one-loop radiative corrections to the coefficients of the NRQCD 
action in a manifestly gauge-covariant manner The coefficient of the er ■ B term in the NRQCD action is 
computed at the one-loop level; the resulting shift of the hyperfine splitting of bottomonium is found to 
bring the lattice predictions in line with experiment. 
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Non-Relativistic QCD (NRQCD) JH is an effective field 
theory that has been applied with considerable success to 
the description of hadrons containing heavy quarks ||2l. 
However, the currently used NRQCD actions do not in- 
clude radiative improvement (with the exception of tad- 
pole improvement), and this will affect the precision with 
which crucial quantities such as the hyperfine splitting be- 
tween the T and the rj\, can be determined [3|. In contrast 
to this, Non-Relativistic QED (NRQED) has been success- 
fully improved and applied to obtain highly precise theo- 
retical predictions for the fine structure of muonium ||4l|5l. 
It is clearly highly desirable to improve NRQCD in a sim- 
ilar manner. This is rendered complicated by the nature of 
the non-abelian gauge interactions in QCD and NRQCD, 
which requires NRQCD to be implemented on a lattice and 
hence makes it necessary to retain the full l/(ma)" de- 
pendences, whereas in continuum NRQED (A/m)" terms 
can be omitted in a consistent manner. Moreover, IR diver- 
gences play a non-trivial role in QCD and NRQCD. 

In this letter, we present the first calculation of radiative 
corrections to coefficients in the lattice NRQCD action us- 
ing the background field (BE) method. We compute the 
one-loop effective action in lattice NRQCD and match it 
term-by-term to the non-relativistic reduction of the one- 
loop effective action for continuum QCD. In particular, we 
determine the one-loop corrections to the coefficient of the 
chromomagnetic cr • B term and the four-fermion spin- spin 
interaction; these corrections are important for the accurate 
calculation of the hyperfine structure of heavy quark states 
using NRQCD. 

THE BACKGROUND FIELD METHOD FOR LATTICE 
NRQCD 

The BE method is a well-established tool to com- 
pute the effective action in quantum field theory. The aux- 
iliary gauge invariance of BEG amplitudes implies that the 
effective action contains only gauge-covariant operators 



which leads to a set of Ward Identities in QCD that reduce 
the amount of calculation necessary to renormalize the the- 
ory. This property is important for operators of dimension 
D > 4 where the loss of gauge-covariance would lead to 
a proliferation of additional operators and is vital to the ra- 
diative improvement of NRQCD which is a non-relativistic 
expansion on operators of increasing dimension; only BEG 
will guarantee the gauge covariance of the improved ef- 
fective action. Whilst the presence of gauge non-covariant 
finite terms with D > 4in the effective action is not per se 
incorrect, they obscure the underlying gauge symmetry and 
greatly complicate the theory and simulation. An attempt 
to match without using BEG would lead to the appearance 
of ultraviolet logarithms, which would have to be cancelled 
by the contributions from additional non-gauge-covariant 
operators. Although BEG does not guarantee that the coef- 
ficients in the effective action are independent of the gauge 
parameter [TOl, in our case we match between theories us- 
ing on-shell quantities and we explicitly find that the coeffi- 
cients are independent of the gauge parameter in both QCD 
and NRQCD. Moreover, the QED-like Ward identities in 
BEG imply that the one-particle irreducible (IPI) vertex 
functions are finite, and that the coupling g is renormal- 
ized only by the contribution from the gluonic self-energy, 
whereas the BE is not renormalized. This is true both in 
QCD and NRQCD, which allows us to match the theories 
by equating two finite quantities. As a consequence of this 
crucial property of BEG, we are free to use different reg- 
ulators in QCD and NRQCD. In particular, we can calcu- 
late the QCD vertex analytically in the continuum using di- 
mensional regularization, or on a fine lattice and taking the 
continuum limit. The latter is particularly convenient for 
checking the gauge-parameter independence of the result, 
since the analytical calculation becomes rather involved for 
general values of the gauge parameter. 

In the following we denote the perturbative expansion 
for a generic parameter wdLSW = ^„^o ^^^"^ck"- 
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FIG. 1. Feynman diagrams to be computed in both QCD and 
NRQCD for matching the cr ■ B term in the NRQCD action. 



MATCHING THE a B TERM 

The effective action for continuum QCD contains the 
following terms involving the fermion fields: 
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which after renormalization of the first term gives 



with 
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which after renormalization becomes 



V„\^R,i,^^,A\ = c,Zl^Zrz''^i,\ 
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We require that the anomalous chromomagnetic moment in 
QCD and NRQCD be equal and find the matching condi- 
tion 



and at tree level and one-loop order we find 
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The NRQCD contribution to €4^^ contains a logarithmic IR 
divergence |^ log(/ia), which combines with the IR loga- 
rithm from the QCD result above to yield an overall loga- 
rithmic contribution — |^ log(Afa). 

Besides the ordinary diagrammatic contributions calcu- 
lated below, we also need to take into account the con- 
tributions from the mean-field improvement U 1— )• U/uq, 

which affect JZ^"^ and 6Z^-^^^l Perturbatively, Uq = 

(2) 

1 — agUQ , and the contributions from inserting this ex- 
pansion into the NRQCD action can be worked out alge- 
braically. The final result for the one-loop correction to C4 
is then 
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j^NR-d) - 5Zf - dZ^'^^^ , (10) 



^^tad,(i) _ 5^tad,(i) _ A log Ma 



where 6Zx denotes a finite diagrammatic contribution. We 
expect the coefficient C4 to be gauge-parameter indepen- 
dent for on-shell quarks, since it is directly related to the 
hyperfine splitting, which is a physical quantity. 



where the leading correction is of order 0{as) and comes 
from 5Z^ alone. After performing the non-relativistic re- 
duction by a Foldy-Wouthuysen-Tani (FWT) transforma- 
tion, we find that the term relevant for the determination of 
the chromomagnetic moment of the quark is 



(4) 



A straightforward analytical calculation of the Feynman di- 
agrams shown in figure [T](a)-(b) gives 
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at the one-loop level, where fj, is the infrared cutoff. 

The effective action for NRQCD contains the spin- 
dependent term 



r.[V^,V/,A] = c4zrV'^^V 



(6) 



THE FOUR-FERMION SPIN-SPIN INTERACTION 

In NRQCD the hyperfine splitting in the bb system 
also receives a contribution from the spin-dependent four- 
fermion operators generated by QQ — t- QQ scattering in 
the colour singlet channel. It is conventional to write these 
contributions using a Fierz transformation ifTl fTTIl 
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where ip and x the quark and anti-quark fields, respec- 
tively, treated as different particle species with correspond- 
ing representations of their spin and colour algebras. The 
spin-independent contributions to di and d2 from QQ scat- 
tering are not included as they do not influence the hyper- 
fine structure. In QCD the two continuum diagrams are 
shown in figuresjSj^a) and[2];b), and in NRQCD all diagrams 
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in figure |2] need to be calculated. The one -loop contribu- 
tions to the renormalization constants for the operators in 
eqn. ( [TT| | take the form, respectively, 

= -\zn. (12) 

The last term in both expressions is the Coulomb singu- 
larity arising from the Coulomb gluon exchange in figure 
[2]; a). For QCD these expressions were verified numerically 
and for both QCD and NRQCD were shown to be gauge- 
parameter independent; there are two independent colour 
trace combinations, each of which is separately gauge inde- 
pendent. In the numerical calculations we used IR subtrac- 
tion functions to analytically remove both IR and Coulomb 
divergences; this greatly improved convergence. For QCD 
we find 



The matching parameters for the term in the NRQCD ac- 
tion, including the two-gluon annihilation contribution to 
di Ifm . are then 

= -3^2-^(2 -2 log 2) , 

d2 = |^-^Af« + ^ log Ma. (14) 

IMPLEMENTATION AND RESULTS 

To perform the calculation in NRQCD, we employ the 
HiPPy and HPSRC packages for automated lattice pertur- 
bation theory llT2l[T3l . which we extended to deal with the 
modifications of the usual Feynman rules engendered by 
the use of BFG lfT4llT5ll . Specifically, in the expansion of 
a gauge link the background fields must always appear 
to the right of the quantum fields and so not all order- 
ings of background and quantum fields can arise. Addi- 
tional contributions to all purely gluonic vertices including 
exactly two quantum gluons arise from the gauge-fixing 
term and additional ghost field vertices are generated which 
have been included but are not needed for the present cal- 
culation. For further implementation details the reader is 
referred to IT61 . 

For the cr ■ B operator matching we compute the dia- 
grams in figures[T](a)-(f) and for the four-fermion operator 
matching we compute the diagrams in figure |2] We use the 
HPSRC library, which includes a parallel implementation 
of VEGAS ifTTl . as well as routines for automatic differen- 
tiation of Feynman diagrams [ 18 |. We carried out a number 
of checks of the calculation. Firstly, we replicate the known 
IR logs correctly. We find that the coefficients of these logs 
are gauge-parameter independent and, since this is not true 
of the contributions from individual diagrams, it provides 



Ma 


1.95 


2.8 


4.0 




-5.164(7) 


-4.913(6) 


-4.739(6) 


5Zm 


1.512(1) 


1.022(3) 


0.723(2) 


SZ^J"^ 


4.387 


4.077 


3.841 


oZjn 


-1.092 


-0.787 


-0.641 


C4 


0.728(7) 


0.799(7) 


0.842(6) 


di 


0.043(7) 


-0.703(14) 


-1.734(25) 


d2 


-0.060(2) 


0.189(5) 


0.532(8) 



TABLE I. Renormalization parameters of the cr ■ B and the four- 
fermion terms defined, respectively, in eqns. and ( |14[ l. 
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FIG. 2. Feynman diagrams to be computed in both QCD and 
NRQCD for matching the four-fermion terms in the NRQCD ac- 
tion. There are two diagrams with the topology of (c). 

a strong check. Second, we check that the non- logarithmic 
part of the result is similarly gauge -parameter independent 
where the individual contributions are not. For matching 
the four-fermion terms it is vital to employ IR subtraction 
functions to remove logarithmic and Coulomb IR singular- 
ities. 

For NRQCD, we used the action from ||2l with stability 
parameter n = 2, and we used the Symanzik improved 
gluon action fT9\, which were also used by the MILC col- 
laboration [20] whose configurations were used in [2J. We 
find 

-^^r^^^^ = - (i + pi^) ^su'o' (15) 

The tadpole contribution to SZ^^'''^^ comes from the mean- 
field improvement of the improved field- strength tensor 
and from the cross-multiplication of the tree-level a ■ B 
term with the tadpole corrections terms in Hq [2J. The 
overall result is 

;>7tad,(l) _ /l3 I 3 3_\ „,(2) 

"^o- ~ 3 AMa 8n{Ma)2 A(Ma)^ J -^0 ■ 

(16) 

(2) 

We chose the Landau mean link to be nj, ' = 0.750 |i2]J. 
Our results are shown in table U Whilst there is no sub- 
stitute for including these radiative corrections in a simu- 
lation, we note that both operators give a contribution to 
the hyperfine splitting that is dominated by a contact term. 
The leading contribution from the gcr ■ B term, already 
included in the simulation, is 0{as) and so the radiative 
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Ma ay((j'*) 


Correction % 
4-fermion a ■ B 


hfs (MeV) 
ref. ^21 


hfs (MeV) 
corrected 


1.95 0.216 
2.8 0.249 
4.0 0.293 


-1.2(1) +31.4(3) 
+11.9(2) +39.8(3) 
+35.7(4) +49.3(3) 


56(2) 
50(2) 
41(2) 


73(3)(5)(6) 
76(3)(6)(5) 
76(3)(7)(4) 



TABLE II. Corrections to the bottomonium hyperfine splitting re- 
sults of 1^1 arising from the radiative improvement of the action. 
In the last column the errors are statistical, O(a^), and relativistic 
corrections. 



correction to C4 and the leading contribution from the four- 
fermion terms in eqn. (Ill both contribute at O(a^) giving 



an 0{a) correction to the measured tree-level contribution. 
A reasonable estimate for the multiplicative correction to 
the tree-level prediction for the hyperfine splitting is then 

av{q*) (2c^l'> -^{d,-dS] , (17) 
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where we chose q* = vr/a. Applying our results to the 
hyperfine splitting of bottomonium, we find the corrections 
given in table [ll]f or the data points of [2|. On all lattices the 
correction is positive and the remaining 0{a^) error in the 
NRQCD predictions of [2 | is reduced to be within errors. 



CONCLUSION 

In this letter, we have applied the BF method to lat- 
tice NRQCD for the first time and have computed the 
one-loop radiative correction to the coefficient, C4, of the 
cr • B operator and the one-loop radiative contribution to 
the coefficients, di and (^2 of the four-fermion contact op- 
erators that affect the hyperfine structure of heavy quark 
mesons. The gauge independence of our calculation was 
explicitly checked by carrying out both relativistic and non- 
relativistic calculations in the lattice theory. This is pos- 
sible because in BFG all calculations are UV finite. Our 
results are summarized in table [I| and in eqns. (fTT]l and 



(14 1. In particular, in eqn. (|11| there is a negative cor- 
rection to C4 due the IR divergences. However, it turns 
out that the constant terms more than cancel this effect and 
the correction to C4 is positive. Whilst there is no substi- 
tute for including these corrections in a simulation, we have 
given an estimate for the correction to the T — rjh hyper- 
fine splitting measured by Gray et al. [2| in table |ll] The 
result is to reduce the lattice spacing dependence to within 
errors and to give an estimate for this hyperfine splitting 
of 73(3) (5) (6)MeV to be compared with the experimental 
value of 69. 3 (2. 8) MeV [22|. The errors shown are statisti- 
cal, O(a^), and due to relativistic corrections, respectively. 
The elimination of 0{aa^) errors and the agreement with 
experiment gives us confidence that the calculations are ro- 
bust. 



The determination of the one-loop radiative corrections 
to the coefficients of the p^, Darwin and spin-orbit terms 
and other four-fermion contact terms, as well as more de- 
tails of the calculations will be presented in a longer paper 
in the near future. 
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